Long-range correlations in the wave functions of chaotic systems 
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We study correlations of the amplitudes of wave functions 
of a chaotic system at large distances. For this purpose, a 
joint distribution function of the amplitudes at two distant 
points in a sample is calculated analytically using the super- 
symmetry technique. The result shows that, although in the 
limit of the orthogonal and unitary symmetry classes the cor- 
relations vanish, they are finite through the entire crossover 
regime and may be reduced only by localization effects. 



During the past decade, a substantial progress has 
been achieved in the studies of quantum properties of 
classically chaotic systems [Q. At present, we have quite 
a complete information about energy spectra of non- 
integrable systems in the limit of a "hard" chaos as well as 
in the crossover regimes to integrability or localization. 
However, the theoretical background for understanding 
the statistics and structure of wave functions of chaotic 
systems is still a growing up field stimulated by the re- 
cent experiments on electron transport in quantum dots 
under the Coulomb blockade conditions Q , as well as by 
the studies of the microwave irradiation in disordered and 
chaotic cavities ||^ . According to the resonant tunneling 
models of the transport in the Coulomb blockade devices 
iQ,^, the height g of the conductance peaks measured in 
the experiments on the quantum dots can be related to 
the local amplitudes (pai'<^i,r) of the wave function of the 
resonant level as 
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-, where broadening Fi^r oc \ipa (ri.r)r 



is due to the electron escape to the right- and left-hand- 
side electrodes Hence, the statistics of the resonant 
tunneling conductances in this kind of measurements is 
determined by the statistics of the electron wave func- 
tions in the dot. Moreover, the statistics and spatial cor- 
relation properties of microwaves in random media have 
become the subject of direct measurements performed 
already in various parametric regimes ||^. 

Below, we demonstrate the existence of long-range cor- 
relations in the wave functions of a chaotic system sub- 
jected to a weak magnetic field, that is, in the crossover 
regime between the orthogonal (with time-reversal sym- 
metry) and unitary (broken time reversal symmetry) en- 



sembles. This conclusion is based on the analysis of the 
joint probability distribution function /2(pi,P2) which 
describes the statistics of local densities of an eigenstate 
in a chaotic cavity taken at two distant points in the 
sample and its comparison with the statistics of a single- 
point distribution /i(pi). The joint distribution func- 
tions /7v(pij ■■■tPn) of local densities of a wave function 
are defined for an arbitrary N as follows 
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where the angular brackets stand for averaging over irreg- 
ularities in the system or over energies e ||7| in a chaotic 
billiard, V is the volume of the system and A is the mean 
level spacing. The sum in Eq. (yj) is taken over all eigen- 
states, ipa and being cigenfunctions and eigenenergies. 
The joint probability distribution function /2 naturally 
appears in the description of fluctuations of the resonant 
tunneling conductance peaks, g, in quantum dots or can 
be directly studied in the microwave experiments jsj . 

Direct computation of the local density distribution 
function /i for disordered and chaotic systems in the lim- 
iting cases of the orthogonal and unitary ensembles ] 
leads to the well known Porter-Thomas distribution ITC ] 
which is merely a Gaussian distribution of wave function 
amplitudes. The Gaussian distribution has been recently 
confirmed by the numerical studies of high-lying chaotic 
eigenstates . It can also be derived semiclassically as- 
suming that classical orbits cover uniformly the energy 
surface in the phase space The detailed analysis of 
the joint distribution function /2 performed in Refs. p^ ] 
for these two universality classes has shown that intrin- 
sic correlations in the wave functions of a chaotic system 
decay already on several wavelengths and the function /2 
can be written at large distances as 



/2 (pi,P2) -> /i (pi) /i iP2) 



(2) 



where the function /i (p) describes the Gaussian distri- 
bution of local amplitudes (/?(r„), p„ = y |(p(r„)|^. In 
particular, the latter result has been used in Refs. 
upon deriving the distribution function of conductance 
peaks in the orthogonal and unitary symmetry ensem- 
bles which were later verified experimentally, Refs. 
Moreover, one can provide a full analytical description 
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of a decay of correlations of local amplitudes of chaotic 
wave functions in each of the symmetry classes |l3| , p^ . 
The law describing the decay of the correlations can be 
derived either using the supersymmetric cr-model tech- 
nique or with Berry's conjecture p2| about the uniform 
distribution of the classical orbits in the phase space , 
and the result for the orthogonal case is in a good agree- 
ment with that of the microwave experiments Q . 

The suppression of correlations within a chaotic wave 
function at large distances (Eq. (^) seems to be so 
natural for systems in which the semiclassical version 
of quantum mechanics deals with large random phases 
that it was not called in question also in the regime of 
a crossover between the orthogonal and unitary symme- 
try classes, where the explicit calculations have not been 
made so far. For example, in Ref. |15| it was argued that 
the presence of a small magnetic field cannot considerably 
change the law describing the decay of the correlations, 
and this assumption was implicit in other works. 

In this Letter we show that due to the Aharonov- 
Bohm effect, the intrinsic correlations implicit into a 
chaotic wave function do not vanish at large distances 
in the entire crossover regime, although they do in the 
orthogonal and unitary limits (corresponding to zero 
and sufficiently high magnetic fields, respectively). In 
weak magnetic fields, the correlations survive at distances 
r = |ri — much longer than the mean free path I (that 

is, f2{Pl,P2) 7^ fl{Pl)fl{P2))- 

Calculations presented below are performed using the 
supersymmetry technique In particular, the use of 
this method enables one to calculate the one-point dis- 
tribution function fi{p) for an arbitrary magnetic field 
(a somewhat different although related quantity was 
calculated in Ref. The function /2 can be com- 

puted in a similar way. Although the calculations are 
not more difficult than those in Ref. j|] they lead to the 
result about the long range correlations which could not 
be extracted from the function /i alone. 

To carry out the computations within the supersymme- 
try scheme, one should write the quantities to be studied 
in terms of the Green functions defined, as usual, as 

Gf -^(n, r^) = Yl ^a(ri)^:(r2)/ [e - ± . (3) 



L[i;]^i / V(r) e-Ha-U{r)-i-fA/2 V(r)dr (5) 



where Hq is the free particle Hamiltonian and U{r) de- 
scribes the impurity potential. The diagonal matrix 
A has components A^^ = — A^^ ~ 1 and is an 8 
-component supcrvector containing both commuting s 
and anticommuting x elements. The reduction of Eq. 
(^ to the integral over the supervectors ip can be done 
by expanding (5-functions in the Green functions and 
G^ and using standard formulae of Gaussian integra- 
tion. Then one averages over the irregularities and, fi- 
nally, uses the Hubbard-Stratonovich transformation to 
decouple the "interaction" term {tpi/j)'^ obtained after this 
averaging by an integration over the supermatrix field Q. 
This gives us a possibility to integrate over the supervec- 
tors ip and to reduce the following computation to that 
of an integral over the supermatrices Q. Integrating over 
the supcrvector ip we take into account only pairing at 
coinciding points. The approximation is justified pro- 
vided the distance between the points ri and r2 is large 
enough (much larger than the wavelength). This means 
that we consider the asymptotic form of the distribution 
function of local amplitudes of the wave function taken 
at essentially distant points. All manipulations are stan- 
dard llq,Ol, and we arrive at 
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where (...)q stands for the integration over super- 
matrices Q with the free energy F [Q] , and 22 = 



-<2), zi = (0,0,e<Se-<^ 
In the presence of a magnetic field, F[Q] is given by 



(0,0,0,0,0,0,e*^^e-<2), zi = (0, 0, e<S e^^'-i , 0, 0, 0^ 
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It is not difficult to show that the distribution function 
/2, Eq. (0), can be represented as follows: 

/2(pi,P2) - - lim 4(/3 t dt f dr3[Gf(r,r)] 

B^l 9P Jo J 
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S [pi +zi:4G^(ri,ri)] 5 [p2 - zF^Gf (r2, r2)] ), (4) 

where Ya = 7/3yt and Yr = lPV{l - t). Next, we 
express the Green functions entering Eq. (^J in terms of 
integrals over supervectors with the Lagrangian L 



where 0o = hc/e is the flux quantum, ag is the sample 
geometry dependent factor, D is the classical diffusion 
coefficient, Ec is the Thouless energy. Eqs. (HJ^) are 
written under an assumption that the supermatrix cr- 
model describing fluctuations of the supermatrix Q is 
zero-dimensional (OD), which is correct for not very large 
Pi,2. The distribution function /2 can also be calculated 
for arbitrary pi_2 in a way that it has been done for the 
function /i in Ref. [ pO[ . This would require going beyond 
the OD (T-model and is out of the scope of this pubhcation. 

Calculation of the integral over Q in Eq. (^ can be 
carried out using the parametrization suggested in Ref. 



2 



p8[ which separates "cooperon" and "difFuson" variables. 
The quantum limit 7^0 considerably simplifies the in- 
tegration over the difFuson variables since the main con- 
tribution comes from that part of the non-compact sector 
of the order parameter Q where \znQzn\ ~ I/7 ^ oo. In 
this limit, the difFuson degrees of freedom can be inte- 
grated out and one ends up with a version of some kind 
of a reduced a-model adapted for describing the proper- 
ties of a single quantum state ]20| and containing only 
the integration over the cooperon variables, 



A 
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In Eq. (||), iV = 1, 2, functions P and An are defined as 



P = Xic + 2 {r],ri* - KcK*) (Ale ~Xc),A„^l + ^- 
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L - V^A2^- 1 (1 + 277,ry:) (1 - 2k,kI) + 4yT^ry,K„ 

and dRc — dr]cd'r]*dK*dKc, with rjc, ijl, Kc, k* being 
Grassmann anticommuting variables. 

Skipping technical details of the algebraic manipula- 
tions with Eq. (H), we are able to represent the joint 
distribution function /jv, N = 1, 2 in the form 



In 
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B{X,x)'[lM{pn,x)da 
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where 



B{X,x) = [{xX^ - 1)^2{X) + $i(X)] 



M{p,x) — y/x exp{—px)lQ (^p\/x^~~xj 



Here, 10(2) is the modified Bessel function, and 



= e'^'dy, $2(^) 



X2 



In order to establish the presence or absence of corre- 
lations between fluctuations of the values pi,2 = |¥'(i"i,2)| 
at distant points in the sample in the crossover regime, 
the joint probability distribution function /2 and its mo- 
ments have to be compared to the distribution function 
/i. Explicit integration over x in Eqs. (^) shows that 
in each of two limits X — > and X — > cx), f2{pi,P2) 
takes a separable form and satisfies Eq. (0). That 

— p ■ /2 

is, /2(pi,P2) — T^i=i,2^y^= in the orthogonal case and 
/2(pi,P2) = exp(— pi — P2) in the unitary one. However, 
the distribution function is not separable at any finite 
magnetic field, which means that even in the limit of a 
large distance between the points ri and r2 the correla- 
tions do not vanish. 



To illustrate the existence of long-range correlation im- 
plicit in each individual wave function in a quantum bil- 
liard subjected to a weak magnetic field, we calculate the 
correlation function K2,2{X) defined as p]| 

/•OO 

K2,2 (X) = / dpidp2p\pl [f2 {Pl,P2) - fl (Pl) fl (P2)] 

Jo 

Although an analytical expression for K2^2 can be found 
from Eqs. (^, we give here only its graphic represen- 
tation. Fig. 1 illustrates how the correlations evolve as 
function of a normalized flux through the sample area. 
As one can see from Fig. 1, the correlations are present 
all over the crossover regime: from the smallest to high- 
est values of the normalized flux cj). These correlation 
are never strong (of the order of one percent), since the 
maximal value of the correlator max{K2) « 0.05 has 
to be compared to the square of the second moment 
P2 — J dpp^ fi{p) which varies from the value of 3 to 
2 in between of the orthogonal and unitary symmetry 
classes. 

The distant correlations in the wave functions can be 
also traced in the joint probability to find simultane- 
ously two zeros of the wave function ip{vi^2) = 0. At 
small fluxes corresponding to X <C 1, [/i(0)]2 w ^X^"^ 
whereas /2(0,0) « f^-^. Since the crossover between 
the orthogonal and unitary ensembles occurs at rela- 
tively small magnetic fluxes through the sample area, 
4> ~ 0o(A/i?c)^^^, the non-separability of f2{pi,P2) as 
a whole is most pronounced when X ~ 1. In Fig. 2, we 
show the difference between /2 and fi x fi a.t X — 1. 

Starting from the pioneering work by Berry | p2[ , it was 
generally believed that in classically ergodic systems the 
local wave function density can be imagined as a result of 
superposition of an infinite number of plane waves with 
random phases and equal momenta. The randomness of 
the phases suggests the Gaussian randomness of the am- 
plitude ifa and vanishing correlations at large distances. 
Our result means that a weak magnetic flux through the 
sample area introduces some correlations in ipa (r) related 
to the fact that the Aharonov-Bohm phases taken by an 
electron moving in a quantum billiard cannot be arbitrar- 
ily large. The lengths of geometrical paths attributed to 
an electron classical trajectory in a semiclassical picture 
of quantum mechanics are limited, since the semiclassics 
breaks down at the time scale longer than the Heisen- 
berg time, tn ~ h/A. Hence, the encircled magnetic 
field fluxes are also limited, so that the classical ergodic- 
ity does not always lead to the complete randomness of 
the wave function phases. Another, better known illus- 
tration of this fact can be given by a disordered system 
with a ring geometry. In the latter case, the flux through 
the ring can exist even if the magnetic field in the wire 
cross-sectional area is negligibly small, which generates 
regular periodic oscillations of the wave function phases. 

The amplitude correlations discussed above persist all 
over the chaotic quantum billiard. One may ask a ques- 
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tion about what happens if the sample size is much larger 
than the mean free path in a disordered system. In fact, 
the result of Eq. (||) is valid as far as the OD cr-model may 
be used, that is, unless localization effects become impor- 
tant. Hence, in ID and 2D weakly disordered samples the 
distance between the observation points must be smaller 
than the localization length, whereas in a 3D metal the 
distance between them is not limited. The correlations 
can be considerably reduced also for large values of am- 
plitudes originated from prelocalized states. In this case, 
one has to go beyond the OD tr-model |^,^ . 

In conclusion, the joint distribution function of ampli- 
tudes of chaotic wave functions is derived for an arbitrary 
magnetic flux. Its form manifests the long-range spatial 
correlations existing in the entire crossover regime from 
orthogonal to unitary symmetry classes. This indicates 
that due to the Aharonov-Bohm effect the phases of the 
wave functions are correlated even if the corresponding 
classical motion is ergodic, which has not been antici- 
pated in previous semiclassical theories. 
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